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ABSTRACT 

The acoustic behavior of a flexible fibrous material was studied 
experimentally. The material consisted of cyllndrlcally shaped fibers arranged 
In a batting with the fibers primarily aligned parallel to the face of the 
batting. This type of material was considered anisotropic, with the acoustic 
propagation constant depending on whether the direction of sound propagation 
was parallel or normal to the fiber arrangement. Normal Incidence sound 
absorption measurements were taken for both fiber orientations over the 
frequency range 140 to 1500 Hz and with bulk densities ranging from 4.6 to 

3 

67 kg/m . When the sound propagated In a direction normal to the fiber 
alignment, the measured sound absorption showed the occurrence of a strong 
resonance, which Increased absorption above that attributed to viscous and 
thermal effects. When the sound propagated In a direction parallel to the 
fiber alignment, Indications of strong resonances In the data were not present. 
From comparing these two sets of data and from considering the material 
structure, the resonance In the data for fibers normal to the direction of 
sound propagation was attributed to fiber motion. An analytical model was 
developed for the acoustic behavior of the material that displayed the same 
fiber motion characteristics shown In the measurements. 

LIST OF SYMBOLS 

c effective propagation velocity In material, Eq. (A20) 

c 

c adiabatic speed of sound In air 
0 



Cp specific heat at constant pressure 
c y specific heat at constant volume 
d fiber diameter (1 . 254x10^ m) 

H porosity, 1 - (p b /p f ) 
h heat transfer coefficient 

i <v - 5 

K heat transfer parameter, Eq. (A14) 
k complex propagation constant 
k' thermal conductivity of air 
L depth of material 
it length of a fiber 
p acoustic pressure 

p' acoustic pressure In air Incident on material 
P q ambient air pressure 

q fiber length with motlon/total fiber length 
R gas constant 

r resistance coefficient of a fiber 
T q ambient air temperature 
t time 

u acoustic particle velocity In the material 
u^ velocity of fiber motion 
u volume velocity per unit area 
u' acoustic particle velocity In air Incident on 
material 

Vj. volume of a fiber 
W characteristic Impedance 
x distance along x-axIs 
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Z acoustic Impedance 

<* n normal absorption coefficient 
Y ratio of specific heats, 1.4 for air 
e acoustic variation In temperature 
v viscosity of air 

v kinematic viscosity of air 

5^ displacement of a fiber 

p acoustic density of air 
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p Q ambient density of air (1.21 kg/m ) 
bulk density of material 
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p^ density of fiber material (1440 kg/m ) 
a viscous loss term, Eq. (A5) 

w radian frequency 

resonant frequency of a fiber 

INTRODUCTION 

Acoustic propagation In porous materials has been studied since the time 

of Rayleigh. The early approach was to assume that the material was rigid and 

that the sound moved through small pores and lost energy through viscosity and 
1 2 

heat conduction. ' Later, the additional complexity of material motion was 
3—6 

added. Since the material frame Itself was flexible, an elastic wave could 

be sustained In the material and Interact with the wave propagating through the 

air. Thus, additional losses of acoustic energy were made through deformations 

of the material frame. Many of these studies have been done for porous 

materials, In general, while others have been based on modeling a specific type 

7 8 

of material. Attenborough ’ has reviewed most of these studies, so no attempt 
will be made here. This study will concentrate on the behavior of one 
particular type of fibrous material. 
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Two characteristics of acoustic propagation In fibrous materials will be 
considered In this report. The first Is the effect of fiber motion on a 
material's acoustic behavior. In this context, the term "fiber motion" means 
that the fiber moves only under the Influence of a passing acoustic wave. The 
ends of the fiber are assumed to be fixed In space at the joints of the complex 

g 

Interconnection of fibers that make up the material. While Attenborough has 

expressed that no real material has shown the behavior described as fiber 
g 

motion, Kawaslma has analytically studied a model where fiber motion was 
assumed to be present In the material. In this report, fiber motion Is used 
to describe some of the behavior seen In the data. When frame motion Is 
discussed, It refers to the motion of the fibers under the Influence of other 
fibers. An elastic, wavellke motion then occurs through the Interconnected 
frame of the material. Host analysis has been done with this flexible frame 
approach, and much data have been collected to verify the analysis. 

The second acoustic characteristic to be considered here results from the 
forming of a fibrous material Into a blanket. When the fibers are manufactured 
Into a blanket, the fibers are primarily aligned In two-dimensional parallel 
planes, with some fibers arranged across these planes to hold the blanket 
together. This general alignment of the fibers creates an acoustic medium that 
Is anisotropic, a characteristic of the medium where the acoustic propagation 
constant depends on the direction of propagation. Inside a fibrous material, 
the acoustic propagation constant depends on whether the sound Is propagating 
parallel to the planes Into which the fibers are arranged or whether the sound 
Is propagating normal to these planes. 

This study considers a fibrous bulk material called Ke»lar. When It Is 
used as a lining for wind tunnel test sections and for aircraft engine Inlets, 
this material can withstand the severe environmental conditions that may occur 
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without breaking down and dispersing In the flow. Kevlar has also been the 
subject of previous work, 10-12 the results of which are used for comparison 
with the results of this study. 

Data consisting of normal Impedance measurements were taken from a series 

of single-density Kevlar samples, all with the same dimensions. The densities 

3 

ranged from 4.6 to 67 kg/m , and data were taken at Identical densities In both 

directions of sound propagation: normal to the fibers and parallel to the 

fibers. After a description of the experimental apparatus and procedure, the 

results for the normal fiber measurements and the parallel fiber measurements 

are shown and compared In terms of the normal absorption coefficients. This 

Is followed by descriptions of the models of Hersh and Walker 10 and of 
g 

Kawaslma, which are compared with data. Finally, the two models are 
combined to obtain an Improved model for the normal fiber measurements. 

EXPERIMENTAL APPARATUS 

An acoustic duct, originally described In Ref. 13, was modified for use as 
14 

an Impedance tube. As shown In Fig. 1, one end was closed off by a steel 
plate, and a 120-W compression driver with an acoustic horn was attached to the 
other end. The rectangular duct had Inside dimensions of 3.81 by 10.16 cm. 

The plane wave cutoff for these dimensions was approximately 1700 Hz, which 
set the upper frequency limit of the measurements. To allow access to the 
Interior, the top of the duct was made of removable plate sections. In one of 
these duct plates, two 0.64-cm condenser microphones were mounted with a 
3.81-cm separation. 

The measurements were based on the two-microphone transfer function 

15 

technique described by Chung and Blaser. With the assumption that the 
microphones measured the sound at a point, the microphone outputs were sent to 
a two-channel fast Fourier transform analyzer, which calculated the transfer 
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function between those two points In the duct. Upon completion of this 
measurement, the two microphone systems were switched In position, and the 
measurement was repeated. This switching procedure resulted In two transfer 
functions, which were used to calculate a geometric average transfer function 
that was corrected for any gain and phase difference between the two microphone 
systems . 

In addition to the measurement technique, considerations were given to 
minimizing bias and random errors. 16 The bias errors were minimized by 
locating the microphones close to the sample and by using a small analysis 
bandwidth. The distance of the microphone farthest from the sample was 
11.18 cm, and the analysis bandwidth was 5 Hz. The random errors were 
minimized by maintaining a high coherence between the acoustic source and the 
microphone signals. This Is difficult to do at very low frequencies and when 
the microphone spacing equals a half wavelength. The coherence yas always 
greater than 0.99 for frequencies of 140 Hz and above. The microphone spacing 
was not a concern since the half -wavelength frequency was about 4.5 kHz, well 
above the frequency range of the measurements. 

The Kevlar samples were cut from a low-density blanket In the form of a 
batting. The fibers were layered and lightly needled to hold the batting 
together. In addition, the batting had been treated with Zepel, a fluid 
repellant. The nominal specifications for the batting were a density of 

3 

6.4 kg/m and a thickness of 2.54 cm. Each piece was cut slightly less than 
the cross-sectional size of the duct to minimize any binding or clamping of the 
material at the duct walls. 

After being weighed, the Individual samples were placed through an opening 
In the top of the duct Into a sample holder, a 4-mesh screen attached to the 
hard walls near the end of the duct, which held the samples In place against 
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the hard end of the duct. With the screen In place, the sample holder was 
3.81 cm high by 10.16 cm wide by 10.16 cm deep. An example of the sample 
holder setup for measurements with the fiber planes normal to the direction of 
sound propagation (referred to as normal fiber measurements) Is shown In 
Fig. 2(a). When measurements were made with the fiber planes parallel to the 
direction of sound propagation (parallel fiber measurements), the samples set 
up for the normal fiber measurements were taken out and rotated 90°, as shown 
In Fig. 2(b). This arrangement allowed both the normal and the parallel fiber 
measurements to be taken with the same samples, the same sample holder volume, 
and the same density. At the lowest density, the samples were lightly expanded 
to fill the volume of the sample holder. Subsequently, to Increase the 
density, more Individual samples were compressed Into the volume of the sample 

3 

holder. Thus, the variation In density ranged from 4.6 to 67 kg/m . 

RESULTS FROM MEASUREMENTS 

The transfer function data were used to calculate the specific Impedances 
and the normal absorption coefficients for all the samples that were measured. 
These results, summarized In terms of the normal absorption coefficient, are 
shown In Fig. 3 for the parallel fiber measurements and In Fig. 4 for the 
normal fiber measurements. Each figure shows how the absorption coefficient 
curves (a n versus frequency) evolve as a function of density, which Is 
labeled on the third axis of the figures. The density Is lowest at the front 
of the figure and Increases to the highest density at the back of the figure. 
The frequency range for all the data begins at 140 Hz and ends at 1500 Hz. 

The absorption coefficient curves for the parallel fiber measurements, 
shown In Fig. 3, are smooth with gradual changes taking place as a function of 
frequency and density. The primary absorption mechanisms are viscous losses 
across the fibers and heat transfer effects between the air and the fibers. 
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In comparison to the parallel fiber data, the normal fiber data, shown In 
Fig. 4, have a resonance present In the data. This resonance Is especially 
apparent In the higher density absorption coefficient curves, where a sharp 
peak Is present In the curves. When the structure of the material Is 
considered, this resonance Is attributed to fiber motion. The viscous drag 
and the pressure gradient across the length of a normal fiber create enough 
force to laterally displace the fiber. In the case of the parallel fibers, 
the viscous drag and pressure gradient forces created along a fiber are not 
enough to cause the fibers to compress or buckle. 

The results for the normal fiber measurements are subject to different 
Interpretations. At first Inspection, It looks as If, as the density 
Increases, the broad peak In the lowest density curve Is evolving Into the 
sharp peak In the highest density curve. It will be shown that this Is not the 
correct Interpretation and that the evolution of the sharp peak, which Is 
attributed to a fiber resonance, also creates the dip In the lower density 
curves at the lower frequency end. The fiber resonance peak occurs at a 
frequency just above this dip, and It Is broader and lower In amplitude at the 
lower densities on this plot; therefore, It Is not clearly visible. As the 
density Increases, the fiber resonance peak becomes higher and sharper, and It 
Is always next to the dip. As shown, the dip and then the peak gradually move 
to a higher frequency as the density Increases. 

The data from the normal and parallel fiber measurements show that the 
material Is anisotropic. The first obvious difference Is the apparent 
presence of fiber motion when the sound propagates normal to the fiber 
direction. The second difference In the two directions Is seen In the overall 
levels of absorption. By Ignoring the normal fiber resonance and by 
considering the broad shape and level of the normal fiber and parallel fiber 
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absorption coefficient curves, a comparison over the range of densities 

measured can be made. At the lowest densities, the parallel fiber curves are 

at the same level or lower than the absorption coefficient curves for the 

normal fiber data. For example, In Fig. 3 the lowest density sample has 

absorption coefficients such that 0.16 <a n < 0.67, and In Fig. 4 for the same 

density, 0.16 <<* n < 0.78. As the density Increases, both sets of data shift 

until, at the highest density, the parallel fiber absorption coefficient curve 

2 17 

Is 3.5 to 14 percent higher than the normal fiber curve. Others * have made 
similar measurements of absorption In materials where the results were 
dependent on the direction of sound propagation through the material. As In 
this study, the materials were made In sheets, with the fibers arranged In 
planes parallel to the face of the sheet. With higher density materials than 
that used here, their results showed that the absorption wa c higher In the 
material when the sound propagated parallel to the plane of the fibers. This 
Is consistent with the results shown In Figs. 3 and 4. 

FIBROUS MATERIAL MODELS 
Rigid Fiber Model 

Hersh and Walker 1 ^ developed a model to describe the acoustic behavior of 
fibrous bulk material. It Is based on empirical relations derived from work 
done on pressure drop and energy transfer across bundles of circular cylinders 
or fibers. Because It assumes that the fibers are rigid, the model does not 
consider frame motion or fiber motion. This model Is useful, however, since It 
describes the general behavior of Kevlar and It Includes both parallel and 
normal fiber terms. 

After the development of the empirical relations, the one-dimensional 

equations for continuity, momentum, and state are given In the following form 

2 

for a differential volume of the material, where ud /v < 1 : 
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dp au 

at " p o ax 


( 1 ) 


M 2I> 
p o at " - ax “ oU 


P = V 


( 2 ) 

(3) 


All the variables are acoustic quantities for air, and their definitions are 
given In the List of Symbols. The viscous loss term a and the effective 
propagation velocity c g are derived In appendix A. The effective propagation 
velocity through the porous material takes Into account the effects of heat 
transfer between the air and the material. 

Combining Eqs. (1) to (3) and assuming that all acoustic quantities are 
proportional to exp(1ut - kx) , we can get the wave equation and subsequently 
the propagation constant. 



(4) 


The characteristic Impedance for the material Is defined as the ratio of 
the pressure wave traveling In one direction to the volume velocity In the same 
direction, W = p/(uH). From Eqs. (2) and (4), the normalized characteristic 
Impedance of the porous material Is found to be 



The final step to determine the surface Impedance of the material with a 
rigid backing requires the following boundary conditions at the surface: 


u ' = Hu 

The first equation has pressure continuous across the boundary, and the second 
equation has mass flux continuous across the boundary. With the additional 
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condition that u = 0 at the rigid boundary, the normalized Impedance for a 
sample with finite depth L Is 


Z _ __W_ 

PC P c 
*0 0 *00 

and the normal absorption coefficient Is 


coth(kL) 


( 7 ) 


1 - 


p c 
*0 0 


- 1 


P C 
H 0 0 


+ 1 


( 8 ) 


A comparison between the parallel fiber data and the original model of 
Hersh and Walker using only the parallel fiber viscous loss term Is shown In 
Fig. 5 for the absorption coefficient, the specific acoustic resistance, and 
the specific acoustic reactance. The figure Includes measured data (line with 
circles), results from the model (solid line), and results from a modification 
of the model (dashed line) to be described later In this section. In this 
sample of the data, the comparison Is quite good throughout most of the 
measured density range. The model specific resistance has larger deviations 
from the data than the reactance at the higher densities. At lower densities, 
both resistance and reactance show their greatest deviation from the data at 
the higher frequencies. This Is especially evident In the normal absorption 
coefficient plot at the lowest density. The absorption coefficient Is 
sensitive to small changes In resistance at small values of resistance. Thus, 
the model underpredicts the absorption. 

Frame flexibility, which Is not considered In any detail In this study, 

Is another factor that appears to affect the results shown In Fig. 5. As 
mentioned In the Introduction, this Is the case where the material frame 

supports an elastic wave which Interacts with the acoustic wave In the air. 

1 8 

Kosten and Janssen have shown for a general porous flexible material that 
frame flexibility Is Important at low frequencies. Their results show that at 
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some low frequency the absorption coefficient for the flexible frame material 
peaks at a higher level than the absorption coefficient for the same material 
with a rigid frame. As the frequency Increases, the absorption coefficients 
for the two frames tend towards agreement. An example of this behavior Is seen 
In Fig. 5(b). The experimental data (line with circles), which Is assumed to 
Include frame flexibility effects, has peak absorption at about 850 Hz. This 
absorption Is higher than that shown for the rigid fiber model represented by 
the solid line. The two curves, then, tend towards agreement as the frequency 
Increases. Similar behavior Is seen In Figs. 5(c) and 5(d). 

Since the model of Hersh and Walker Is based on some empirical 
correlations, an attempt was made to see If any Improvements could be made In 
the model by varying the fit parameters within the correlation, specifically 
the parallel fiber viscous loss term originally given by Hersh and Walker as 

f p (l - H) =» 1.0 £ 3 . 94(1 - H)°* 413 [l + 27(1 - H) 3 ]j (9) 

The term f p (l - H) was used as a fit parameter to adjust the model Impedance 

to best fit the data Impedance for each density. The results of this fit are 

shown In Fig. 6. The fit for data below 1 - H = 0.02, or a bulk density of 
3 

about 29 kg/ro , Is scattered, and no fit equation comparable to Eq. (9) Is 
possible. However, the fit data above this point do provide a good correlation 
represented by the equation 

fp(l - H) = 1 .8 ^3. 94(1 - H)°' 5 [l + 27(1 - H) 3 ]j (10) 

When this fit Is used In the acoustic model, the results are as shown In 
Fig. 5. The lack of significant change between the original parallel fiber 
model (Eq. (9)) and the modified model (Eq. (10)) when compared with the data 
Is most likely due to the lack of Inclusion of frame flexibility effects. This 
Is especially evident at the lower densities where, as shown In Fig. 6, the 
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fit of fp(l - H) did not agree with the rigid fiber correlations given by 
Eq. (9) or Eq. (10). 

Comparison between the normal fiber model of Hersh and Walker and the 
data Is shown In Fig. 7 for the absorption coefficient, the specific acoustic 
resistance, and the specific acoustic reactance. The trends In the data are 
similar to the model; however, the model generally underpredicts the measured 
Impedance. Also, the model shows no Indications of the resonances seen In the 
three higher densities In Fig. 7. The correlation for the normal fiber viscous 
loss Is given by 

f n (l - H) = 0.44 ^16(1 - H)°' 5 [l + 14.75(1 - H) 3 ]j (11) 

Its similarity to Eqs. (9) and (10) shows that the expected behavior of the 
normal fiber model Is the same as the behavior shown for the parallel fiber 
model. The difference would be due to the different levels of viscous loss. 

It Is Interesting to note that for 1 - H < 0.05, which Is the case for all the 

3 

measured densities, (1 - H) « 1 with the result that Eq. (11) Is 

approximately equal to Eq. (10). 

11 12 

The works of Lambert and Smith and Parrott have shown data on the 
Impedance of Kevlar samples that have the same characteristics as shown In 
Fig. 7. Lambert's data for a porosity of 0.94 and a sample depth of 12 cm 
showed a resonance at about 950 Hz. He states that this may be evidence that 
the sample Is vibrating and absorbing energy. Smith and Parrott took Impedance 
measurements for Kevlar with and without a Zepel fluid repellant treatment. 
Their data for 1-cm-thlck Kevlar samples at a porosity of 0.94 showed a 
resonance between 700 to 800 Hz. No other resonance was seen In the data out 
to 3.5 kHz. They attributed the resonance to the possibility of fiber motion. 
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Fiber Motion Model 


g 

Kawaslma developed a model that Included the effects of fiber motion. 

He assumed that each fiber could be represented as a string fixed at both ends. 
To simplify the analysis, the string of length 8. was modeled as a rigid bar 
with the center part of length qt vibrating linearly and with the end parts 
of total length q(l - l) fixed. The vibrating force was the oscillating 
viscous fluid, and the restoring force was governed by Hooke's law. From a 
dynamical analysis comparing the vibrating string to the vibrating bar, q was 
found to be equal to 0.811. The following equations were developed: 


0£ 

at 


9U , 1 - H 3u f 

p o ax ~ p o q H ax 


au 
p o at 




3U f 2 
p f at * p f u f*f = 


dX 


- ^ ( u f - u) 


( 12 ) 

(13) 

(14) 


Equation (12) Is the continuity equation. The change In mass In a 
differential volume Is governed by both the movement of air and the movement of 
fibers. The latter Is represented by the second term on the right side of the 
equal sign. This equation Is derived In appendix B. 

The momentum equation Is shown In Eq. (13). It shows that the motion of 
an elemental mass of air Is driven by the pressure gradient, with damping 
across both the fiber parts that move and those that are fixed. The viscous 
loss coefficient, given by the ratio r/v^. Is assumed to be the same for 
both moving and fixed portions of the fibers. 

The third equation governs the fiber motion. The left side of Eq. (14) 

Is a mass-spring equation for the fiber In the elemental volume, with the 
restoring force governed by the resonant frequency u^. The forcing terms. 
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shown on the right side of the equation, are due to the pressure gradient and 
the viscous drag on the fiber. 

When the acoustic pressure and density are related by the adiabatic speed 

2 

of sound In air (p = c Q p) , Eqs. (12) to (14) can be combined to derive a wave 
equation for the pressure. Thus, the propagation constant Is determined to be 



The characteristic Impedance for the material In this case Is defined by 
using the following equation for the volume velocity per unit area: 

0 = Hu + q(l - H)u f (16) 


This equation shows that the volume velocity Is affected by those fibers that 
move. The result Is that, depending on the phase relationship between u and 
u^, the volume velocity Is higher or lower than that which would be found In a 
rigid material having the same porslty. By using Eq. (16), the normalized 
characteristic Impedance Is found to be 


W 


p o c o 


u 


p o c o 



(17) 
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Using the boundary conditions given In Eq. (6) and the additional fact 
that u • Uf i 0 at the rigid boundary yields the specific acoustic Impedance 
at the surface: 


p o c o p o 


1 + q 


1 - H 
H 


Nil 

n * 


J 


coth(kL) 


where 


W_ 
c o . 


[ H p 0 “ J p o“ 


( 18 ) 


"r/Vf 

_p 0 “ 


+ 1 


fl 


o J 


♦ Cl 


1 - H r/V f 


H 


P 0 “ 


It should be noted that the ratio N'/D' equals the ratio between the fiber 

velocity and the air velocity u f /u. The absorption coefficient Is then 

calculated from the Impedance by using Eq. (8). This model development Is 

19 

similar to that given by Ingard, except that material frame motion was used 
there Instead of fiber motion. 

In Kawaslma's model, the viscous loss coefficient r/v^. Is computed for 
Stokes flow around a long ellipsoid of gyration. The term Is a function of the 
length and diameter of the fiber, and not a function of frequency. No attempt 
was made In this study to calculate this term; It was used as a parameter to 
fit the model to measured data. 

A comparison between Kawaslma's model and the data Is shown In Fig. 8. 

The resistance, reactance, and absorption coefficient curves for both model 
and data follow the same behavioral trends. The obvious feature Is the 
resonance due to fiber motion. The assumption that viscous loss Is not a 
function of frequency Is shown to be Inaccurate, since the viscous loss 
coefficient must provide the loss at each frequency and be a damping term at 
the resonance frequency. The model provides Insufficient damping at the 
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resonance, and It Is much less accurate In predicting Impedance at the low 
frequencies. 

Combined Model 

Two factors that contribute to the Inaccuracy of Kawaslma's model are the 
use of the non-frequency-dependent viscous loss coefficient r/v^ and the 
use of the adiabatic speed of sound for the air Inside the fibrous material. 

It Is proposed here to modify Kawaslma's model by replacing these two factors 
by the equivalent factors from the model of Hersh and Walker. 

The first factor to replace Is the viscous loss coefficient r/v^. If 
there Is no fiber motion, then u^ = 0. Substituting this value Into Eq. (13) 
and reducing results In 


3u d£ 1 - H r_ 

p o at 3 " ax H Vf 


(19) 


Comparing this momentum equation to Eq. (2) shows that for the two equations to 
be equivalent 



Thus, the constant r/v f Is replaced by a function that Is dependent on both 
frequency and porosity. 

The adiabatic speed of sound Is the other factor to replace In Kawaslma's 

model. This Is done by using Eq. (3) to define the relationship between 

pressure and density. When the propagation constant and characteristic 

Impedance are determined, the result of this change Is to replace c by c 

o e 

In Eq. (15) and to multiply the right side of Eq. (17) by (c /c ) 2 . Thus, the 

c u 

effects of heat transfer between the air and the fibers have been taken Into 
account. 

The final calculations, then, for the specific Impedance are done by using 
Eq. (18) with the proper substitutions: Eq. (20) Is used to replace the 
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viscous loss coefficient, and the modified Eqs. (15) and (17) are used for the 
propagation constant and the normalized characteristic Impedance, respectively. 

The results of calculations using this combined model are shown In Fig. 8. 

The only factor unknown a priori In these calculations was the fiber resonant 

frequency All other factors were based on the original models of Hersh 

10 9 

and Walker and Kawaslma. As shown In the figure, the calculated specific 
reactance agrees quite well with the data, but the model does underpredict the 
specific resistance. Factors such as losses due to Internal damping of the 
fibers or due to fibers rubbing together and frame flexibility were not 
considered. This may have led to the Inaccuracies In the model at the lower 
densities, at the lower frequencies, and at the resonant frequency. 

CONCLUDING REMARKS 

Experiments were conducted to study the acoustic behavior of a fibrous 
material. The fibers in this material were arranged such that the sound could 
propagate In a direction normal to the fiber arrangement or parallel to the 
fiber arrangement. This permitted an evaluation of the anisotropic acoustic 
characteristics of the material. 

The results of the measurements were summarized In three-dimensional plots 
of absorption coefficient versus frequency as a function of the sample density. 
When the sound was propagating normal to the fibers, the measurements Indicated 
a resonant condition, which was attributed to fiber motion. This resonance was 
shown to change in frequency as the density Increased. When the sound was 
propagating parallel to the fibers, the measurements showed no Indications of a 
fiber motion resonance. Other than the fiber motion effects, anisotropy was 
shown to exist over the rest of the frequency range of the measurements. This 
was shown by the level of absorption that took place. At low density, the 
parallel fiber absorption coefficients were less than the normal fiber 
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absorption coefficients. The opposite was true at the higher densities. 
There, the normal fiber absorption coefficients were less than the parallel 
fiber absorption coefficients. 

Models were used to confirm these results. The parallel fiber 
measurements were compared with the parallel fiber model of Hersh and Walker. 
Though an attempt was made to Improve the model by adjusting constants In the 
model, no significant Improvement was obtained. For the normal fiber model. 
Improvements were made by combining the model of Hersh and Walker with the 
model of Kawaslma. This combined model accounted for the effects of fiber 
motion, and It accounted for the frequency and porosity dependence of the 
viscous loss term and of the heat transfer parameter. Results for this model 
were found to agree with the data except at the lower densities, where frame 
flexibility, which was not Included In any of the models, had a large effect. 



APPENOIX A - DERIVATION OF EMPIRICAL PARAMETERS 
Viscous Loss Term a 


The one-dimensional pressure drop per unit distance across a bundle of 
fibers numbering N Is 

Afi N ~ 


(Al) 


AX ■ AXA u 

where D Is the fluctuating drag of each fiber and A Is the cross-sectional 
area of the duct. The minus sign assures that the pressure decreases across 
the bundle of fibers. The fluctuating drag Is given by the product of the 
cylindrical area A^ of a fiber and the area-averaged fluctuating shear stress 


AV 


where t av Is approximated by 


D = Vav 


(A2) 


t AV 



(A3) 


The function f Isa nondlmenslonal term that accounts for the average drag of 
a circular cylinder In a bundle of cylinders. It Is a function of porosity. 
Using the definition of porosity as the ratio of the volume of air to the total 
volume of the material yields the following equation: 


1 - H = 


AXA 


Combining Eqs. (A2) to (A4) Into Eq. (Al) yields 

1(1 - H)fu = - au 


Afi 

AX 


- ( 
Id 2 / 


(A4) 


(A5) 


where the following relation was used for the ratio between the circular area 
of a cylinder and Its volume: 

A f 4 

f 4 (A6) 


v f - d 


The viscous loss term a has been defined In Eq. (A5). It Is made 
nondlmenslonal by dividing by p 0 and «. Thus, 
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*>o" 


P 0 <rt 


2 (1 - H)f - 4 


(i) 


(1 - H)f 


The function f Is defined from correlations derived by Hersh and Walker 
For parallel fibers, f 
Eq. (11). 


f In Eq. (9); and for normal fibers, f * f R In 


(A7) 

10 


Effective Propagation Velocity c g 

For a unit volume of material AxA, the rate of Increase In energy per 
unit volume of air minus the net rate of energy out of the volume through the 
air must equal the energy flow from the fibers to the air per unit volume of 
air. This statement can be written as 


_ 30 fo 0£ 

p o v at " p 0 at 

Combining an equation for the heat flow 

Q = hNA f e 

and the definition of the Nusselt number 


e 


(A8) 


(A9) 


Nu 


hd 


(A10) 


we get 


Q = 


Nu k'NA f © 


(All) 


To put Eq. (All) on a basis of per unit volume of air, we divide the equation 
by AxAH. Then, by using Eqs. (A4) and (A6), an equation for e Is obtained: 


e = 


AxAH 


4k 


'Nu ( } - H \ 

2 V H / 


e 


(A12) 


Thus, with this equation, we can rearrange Eq. (A8) to get 


where 


— + — A - f v 11 — ^ 

« *„ 9 ■ <T " ') o 0 9t 


P n RT 

— — = — = (Y - 1)T 
Vv c v 0 


(A13) 
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and 


4k 

d 2 c 


H't") 


Nu 


(A14) 


Equation (A14) can be nondlmenslonallzed as follows: 


c , 

-B 


1 /I - H 


p 0 “ 


, C v/V wC pA p o/<od 2 


H 


Nu 


il /_« 

‘ Pr U 2 


:XV) 


Nu (A15) 

where Pr Is the Prandtl number for air. The correlation for Nu given by 


Hersh and Walker 10 is 

Nu = 1.0 ^5.4(1 - H) 0 5 [l + 3.94(1 - H) 3 ]j (A16) 

Equation (A13) may now be solved for 0 by assuming harmonic motion 
exp(lut) . 


e = 


My-DT, 



(A17) 


The equation of state Is given In Its Instantaneous form as 

(P 0 + P) = (p 0 + p) R(T 0 + e) (A18) 

When only linear terms are retained, this equation becomes 

p = p Q R0 + RT Q p ( AT 9 ) 

and upon substituting In Eq. (A17), a relationship between p and p Is 

obtained. 


P 



where 


C o = RT rt 

e o 


K 

V + 

lY 

K 

1 


. p o w * 

_ 


(A20) 


Note that when K/p Q « » 1 (very low frequencies), c‘ = RT Q , the Isothermal 

2 

speed of sound. Also, when K/p « « 1 (very high frequencies), c = yRT , the 

0 c o 

the adiabatic speed of sound. 
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APPENDIX B - DERIVATION OF EQUATION (12) 

Equation (12) Is derived from two continuity equations. The first 
continuity equation Is for the air In the material, and the second continuity 
equation Is for the fiber material. 


d . ... u au 

it (p o + p )H = “ p o H ax 

3 _ 

at 


P f ( 1 -• H) = - p f q(l - H) 


dU| 

ax" 


Expanding these two equations and retaining only linear terms, we get 

au 

'o" ax 


Iti . u ** 
at 


^o at + H s P « H av 


!ff. 3H H ^ an HI ^ 

at ' p f at ' H at a " p f q(1 ~ H * ax 


(Bl) 

(B2) 

(B3) 

(B4) 


By assuming that the compressibility of the fiber material *s very small 
compared to air and may be neglected, then ap^/at = 0 and Eq. (B4) becomes 

it ■ 1,(1 - H) 5T 

Substituting Eq. (B5) Into Eq. (B3) and rearranging, we get Eq. (12). 


(B5) 


dp au 

— r- = - o - — 


at 


o ax 


p o q 


1 - H 3U f 
H ax 


(B6) 
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FIGURE 1. - ACOUSTIC IMPEDANCE MEASUREMENT DUCT. INSIDE DIMENSIONS. 3.81 BY 
10.16 cm. 
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(B) FIBER DIRECTION PARALLEL TO DIRECTION OF SOUND PROPAGATION. 


FIGURE 2. - MATERIAL SAMPLES PLACED IN DUCT. 
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